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The angular distribution of upward-going muons produced by atmospheric neutrinos in the rock below the 
MACRO detector show anomalies in good agreement with two flavor —> Vr oscillations with maximum mixing 
and Am'^ around 0.0024 eV^ . Exploiting the dependence of magnitude of the matter effect on oscillation channel, 
and using a set of 809 upward-going muons observed in MACRO, we show that the two flavor //^ Va oscillation 
is disfavored with 99% C.L. with respect to Vy, — > z/r- 
PACS 14.60Pq; 14.60Lm 



1. INTRODUCTION 

Neutrino oscillations were first suggested by 
B. Pontecorvo in 1957 after the discovery of the 
if" ^ if transitions. Subsequently, evidence 
for the existance of neutrino oscillation in na- 
ture has been provided by the Superkamiokande, 
Soudan2 and MACRO experiments, each of which 
has presented data which strongly favor atmo- 
spheric neutrino oscillations, in the form of 
disappearance 

The two neutrino oscillation probability in vac- 
uum is given by: 
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where Am^ = - m|( eV^), L(km), £;(GeV), 6 
is the mixing angle and L is the path length be- 
tween the neutrino production point and the lo- 
cation at which the neutrino flavor is measured. 
This simple relation should be modified when a 
neutrino propagates through matter and when 
there is a difference in the interactions of the two 
neutrino flavors with matter ||]. The neutrino 
weak potential in matter is: 
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where the upper sign refers to neutrinos, the lower 
sign to antineutrinos, Gp is the Fermi constant, 
ns the baryon density, Yn the neutron and Vg 
the electron number per baryon (both about 1/2 
in common matter). The weak potential in mat- 
ter produces a phase shift that will modify the 
neutrino oscillation probability if the oscillating 
neutrinos have different interactions with matter. 
Therefore, the matter effect could help to discrim- 
inate between different neutrino oscillation chan- 
nels. According to equation, ( ||) matter effects 
in the Earth could be important for i/^ and 
for the — > Vs oscillations, while for — > Vr 
oscillations there is no matter effect. For particu- 
lar values of the oscillation parameters the matter 
effect increases the oscillation probability, leading 
to 'resonances' (e.g., the MSW effect). 

— > Vg oscillations have been suggested Q to 
explain some features of the atmospheric neutrino 
anomaly. Under most current models, a fourth 
(sterile) neutrino is necessary to explain all the 
reported neutrino anomalies (solar, atmospheric 
and LSND [Q ). Matter effects are important 
§ when £;^/|Am2| > 10^ GeV/eV^, therefore in 
particular for high energy events. The primary 
purpose of this letter is to compare the MACRO 
high energy neutrino events sample with the pre- 
dictions, considering matter effects in the case of 
v^l, Vs oscillations. In MACRO, neutrino oscil- 
lation is observed in three different event topolo- 
gies, having different characteristic ranges of par- 
ent neutrino energies. So-called Up Through 
events are associated with muons which pen- 
etrate the entire detector. The parent neutrinos 
in these events have a median neutrino energy 
around 50 GeV. Internal Up events and Internal 
Down events, together with Up Stop events, Q 
are associated with muons having a track termi- 
nus located within the MACRO detector. The 
parent neutrinos in these events have a signifi- 
cantly lower median energy, of around 4 GeV. 
In this paper, we focus on the high energy {Up 
Through) data sample. A similar analysis has 
been recently published by the Superkamiokande 
collaboration IHI. 



2. DATA ANALYSIS 

The MACRO detector [§. is located in the 
Hall B of the Gran Sasso Laboratory, with a 
minimum rock overburden of 3150 hg/cm^. It 
is in the general form of a large rectangular box, 
76.6 m X 12 m X 9.3 m, divided longitudinally 
into six supermodules, and vertically into a lower 
part (4.8 m high) and an upper part (4.5 m 
high). The active detection elements are planes of 
streamer tubes for tracking and of liquid scintil- 
lators for fast timing. The lower half of the detec- 
tor is filled with trays of crushed rock absorbers 
alternating with streamer tube planes, while the 
upper part is open and contains electronics racks 
and work areas. 
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Figure 1. Zenith distribution of the flux of upgo- 
ing muons with energy greater than 1 GeV for the 
combined MACRO data. The shaded region shows 
the expectation for no oscillations with the 17% nor- 
malization uncertainty. The lower line shows the pre- 
diction for an oscillated flux with sin 29 — 1 and 
Am^ = 0.0024 eV^ 



The Up Through muon tracks we focus on 
in this study come from interactions in the 



rock below MACRO. In these events, the muon 
crosses the entire detector requiring that E^> 1 
GeV. The time information provided by the scin- 
tillator counters determines the flight direction 
of the muon, allowing Up Through events to 
be distinguished from the much more common 
down-going muons. The measured muon velocity 
is calculated with the convention that down go- 
ing muons have /3=velocity/c=-|-l while up going 
muons have (3 = —\. In the Up Through event 
sample, almost 50% of the tracks intercept 3 scin- 
tillators planes. In this case, there is redundancy 
in the time measurement, and (3 is calculated from 
a linear fit of the times as a function of the path 
length. Tracks with a poor fit are rejected. Up- 
ward going muons are selected by requiring that 
the measured velocity lie in the range — 1 .25 < 1//3 
< -0.75. 

The data used in this study have been collected 
in three periods, with different detector configu- 
rations, starting in 1989. The statistics is largely 
dominated by the full detector run, started in 
May 1994 and ended in December 2000 (live time 
5.51 years). The total live time, normalized to 
the full detector configuration, is 6.17 years. 

Several cuts are imposed on the data to remove 
backgrounds caused by radioactivity or showering 
events which may result in bad time reconstruc- 
tion. The primary data selection in this regard re- 
quires that the position of a muon hit in each scin- 
tillator, as determined from the timing within the 
scintillator counter, agrees within ±70 cm with 
the position indicated by the streamer tube track. 
This eliminates events with significant errors in 
timing. In addition, downgoing muons which 
pass near or through MACRO may produce low- 
energy, upgoing particles, which could appear to 
be neutrino-induced upward throughgoing muons 
if the down-going muon misses the detector [|lO| . 
In order to reduce this background, we impose a 
cut requiring that each upgoing muon must cross 
at least 200 g/cm^ of material in the bottom half 
of the detector. Finally, a large number of nearly 
horizontal {cosd > —0.1), but upgoing muons 
have been observed coming from azimuth angles 
(in local coordinates) from -30° to 120°. In this 
direction the overburden is insufficient to remove 
nearly horizontal, downgoing muons which have 
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Figure 2. The MACRO 90% confidence level regions 
computed using the angular distribution only (dashed 
line) and the angular distribution combined with the 
normalization (continous line). 



scattered in the mountain and appear as upgoing. 
We exclude this region from our data. 

After applying the data selections described in 
the previous paragraph, we observe 863 events 
with measured velocities in the range —1.25 < 
< —0.75. Based on events outside the upgo- 
ing muon peak, we estimate that there are 22.5 
background events in this data sample. In addi- 
tion , we estimate that there are 14.2 events which 
result from upgoing charged particles produced 
by downgoing muons in the rock near MACRO. 
Finally, it is estimated that 17 events are the re- 
sult of interactions of neutrinos in the bottom 
layer of MACRO scintillators. After subtracting 
these backgrounds to the UpThrough data set, 
the number of upgoing throughgoing muons inte- 
grated over all zenith angles is 809. 

In the simulation of our upgoing muon data, 
we have used the neutrino flux computed by the 
Bartol group and the GRV94 ^ parton dis- 
tribution set, which increases the upgoing muon 
flux by -1-1% with respect to the Si |l^ parton 
distribution that we have used in the past. For 
low energy channels (quasi-elastic and 1 pion pro- 
duction) we have used the cross section in [l4| . 
The propagation of muons to the detector has 



Figure 3. Reduction factor for sin^26 — 1, two values 
of Am^ and i/„ — > i^s or Un — » Ut 



been done using the energy loss calculation by 
Lohmann et al. ||l^ for standard rock. The to- 
tal systematic uncertainty in the predicted flux of 
upgoing muons, adding in quadrature the errors 
from the Bartol neutrino flux, the neutrino cross- 
section, and muon propagation, is estimated to be 
±17%. This theoretical error in the predicted flux 
is mainly a scale error that does not change the 
shape of the angular distribution. Assuming no 
oscillations, the number of expected events inte- 
grated over all zenith angles is 1122, giving a ratio 
of the observed number of events to the expec- 
tation of 0.72 ±0.026(stat) ±0.043(systematic) 
±0.12(theoretical). 

Figure ^ shows the zenith angle distribution of 
the measured flux of upgoing muons with energy 
greater than 1 GeV for our full upgoing data sam- 
ple, compared to the Monte Carlo expectation for 
no oscillations, and with a z^^ — > i't- oscillated flux 
with maximum mixing and Am? = 0.0025 eV^. 
The shape of the angular distribution has been 
tested with the hypothesis of no oscillations, nor- 
malizing the total predicted flux to that observed. 
The is 25.9 for 9 degrees of freedom (P=0.2%). 
Under the hypothesis of t'^ — > i^r oscillation, the 
best is 7.1 and is outside the physical region. 
The best in the physical region of the oscilla- 
tion parameters is 9.7 (P=37%) for Am^ of 0.0025 
eV'^ and maximum mixing. Combining informa- 



tion from the angular distribution and the total 
number of events according to the procedure de- 
scribed in [|6|,we obtain a peak probability of 
66% for oscillations with Am^ of 0.0024 eV"^ and 
maximum mixing, while the probability for no os- 
cillations is 0.2%. 

The 90% confidence level regions of the 
MACRO upgoing events are shown in Figure |[ 
The limits are computed using the Feldman- 
Cousins procedure . Figure || shows the results 
obtained using the angular distribution alone, 
and the angular distribution together with the in- 
formation due to the overall normalization. The 
90% confidence level regions are smaller than the 
regions obtained by SuperKamiokande p8[ and 
Kamiokande p9| ] for the upgoing muon events. 
This can be accounted for through the follow- 
ing effects: the different energy threshold (Su- 
perkamiokande has an average energy threshold 
of about 7 GeV, MACRO has 1.5 GeV), the use 
of the Feldman-Cousins procedure, and the fact 
that our best point is outside the physical region. 

3. TWO FLAVORS STERILE NEU- 
TRINO OSCILLATIONS AND TAU 
NEUTRINO OSCILLATIONS 

In the Vfi — Vs oscillation scenario, the matter ef- 
fect changes the shape of the angular distribution 
and the total number of events with respect to 
vacuum oscillations. Large matter effects are ex- 
pected for neutrinos near vertical incidence, due 
to the large neutrino path length in this case, and 
to the increase in the density of the Earth near 
its core. Assuming maximal mixing, as suggested 
by all available data, the matter effect produces 
a reduction of the oscillation effect, and results in 
an upgoing muon flux closer to that predicted by 
the no oscillation scenario. This effect would be 
most pronounced for directions near the vertical 
|6| . Figure || shows the reduction with respect to 
no oscillations for maximal mixing for — s- Vs 
and Vt oscillations, with Arn^ = O.OOley^ 

and Arn^ = O.OleF^. We have tested the shape 
of the observed upgoing muon angular distribu- 



2.5 



sterile neutrino 
neutrino 




tion against the hypothesis of z^^ 
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with maximum mixing. The best is 20.1 with 
9 degrees of freedom. Combining the information 



Figure 4. The ratio between the data in two bins 
(dashed line) and the comparison with the Vs and v-r 
oscillations with sin^26 — 1. The error bar includes 
statistical and systematical error. 



obtained from the angular distribution and the 
normalization the highest probability obtained is 
8% for maximum mixing and Am^ = 0.006eV^^. 
A statistically more powerful test is based on the 
ratio between the number of events in the two 
angular regions cos{9) > 0.7 and cos{9) < 0.4 
as shown in Figure ^ This quantity is statisti- 
cally more powerful than the in 10 bins be- 
cause data are binned to maximize the difference 
between the two hypotheses to be tested and be- 
cause the ratio is sensitive to the sign of the vari- 
ation (while the is not). In addition, this ratio 
is insensitive to most of the errors in the theo- 
retical prediction of the v flux and cross section. 
The primary disadvantage of this statistic is the 
loss of some features of the angular distribution. 
We have chosen slightly different angular regions 
than suggested in the original proposal for this 
statistic, presented in ref [Q. In our study, the 
angular regions used are based on a Montecarlo 
study of the intervals providing the maximum dis- 
crimination between the v^, — > Vs and the i/^ —^ Vt 
oscillation hypotheses. 

The ratio of the flux of upgoing muons in two 
angular intervals is insensitive to uncertainties in 
the overall i' flux and cross section, as pointed out 



in the last paragraph. Several effects do, how- 
ever, lead to systematic errors in this ratio. For 
example, uncertainties in the K^tt fraction in at- 
mospheric air showers, and the different angular 
distributions of neutrinos produced by these par- 
ents, leads to approximately a 3% systematic er- 
ror]^ in the predicted value for this ratio. An- 
other theoretical error, at the level of approxi- 
mately 2% for MACRO, results from uncertain- 
ties in the neutrino cross sections, and the dif- 
ferent energy distributions of neutrinos arriving 
from the horizontal and vertical directions. A fi- 
nal source of systematic error in the prediction 
of the flux ratio results from the seasonal varia- 
tions of the atmosphere's density profile, and the 
fact that the neutrino flux is computed for the 
standard United States atmosphere jll] not tak- 
ing into account variations of the density profile 
with latitude. Seasonal variation of the high en- 
ergy muon flux has been observed by MACRO 
at 42° North latitude, where a 3% difference was 
observed between summer and winter. At more 
extreme latitudes, Amanda l23|, which operates 
near the South Pole, observes a 20% difference 
between winter and summer. A precise estimate 
of the seasonal variation of the high energy neu- 
trino flux is rather difficult to obtain because it 
requires knowing the density profile of the atmo- 
sphere over the entire Earth. We have performed 
a simplified estimate of the size of this effect based 
on an analytic neutrino flux calculation [ p4[ and 
the CIRA-86 atmosphere tables According 
to this calculation the amplitude of the seasonal 
variations of the ratio of the vertical to horizontal 
neutrino flux is of the order of ±2.6%. Assuming 
a sinusoidal variation during the year, this ampli- 
tude corresponds to a root mean square value of 
about 1.3%. Dividing the MACRO data into a 
winter set (including the months from November 
up to April) and a summer set (the remaining 
months), we observe a difference in the ratio of 
the flux in the two angular bins of 19% ± 17% 
between the two data sets, with a smaller value 
in the summer as expected for the seasonal vari- 
ation, compatible inside the large errors with the 
expectations. We include in our estimate of the 
total systematical error in the predicted flux ra- 
tio a 1.3% contribution due to seasonal variations. 



The systematic error due to the fact that the neu- 
trino flux is calculated using the standard United 
States atmosphere has been estimated to be less 
than 1%. Accounting for all contributions to the 
systematic error, we estimate that the total un- 
certainty in the predicted value for the flux ratio 
is 4%. 

The total systematic error in the measured 
value of the flux ratio has been estimated to be 
4.6%. This error is due to uncertainties in the 
efficiency of the analysis cuts and detector effi- 
ciencies; it could be reduced in the future with a 
reprocessing of the data to correct for the change 
of the apparatus operating conditions with time. 
Combining in quadrature the theoretical error 
and the experimental error we obtain a total error 
in the ratio of about 6%. 

In the full upgoing muon data set, there are 305 
events with cos{9) < —0.7, and 206 events with 
cos{9) > —0.4, giving a value for the flux ratio 
of Re^p = 1.48 ± 0.13,tat ± 0.10,2,,*. This mea- 
sured value can be compared with 
and Rsterilernin='2-^&, which are the minimum 
possible values of R for — > v-r and — os- 
cillations respectively, for maximum mixing and 
Am^ of 0.0025 eV'^. For values of sin'^2e < 1 the 
value of R is larger than Rmin both for i/^ — z/, 
and — > v-r- We note that this ratio does not 
have a gaussian distribution - the errors are re- 
ported only to give a crude estimate of the sta- 
tistical significance. The corresponding one sided 
probability PbestT of measuring a value smaller 
than Rexp, assuming a true value for the ratio 
of RTmim is 8.4%. For — Vg the probability 
Pbestster is 0.033%. The ratio of the probabili- 
ties PhestT I Pbestster is 254. This implies that 
i^fi.—Vs oscillation (with any mixing) is excluded at 
about 99% C.L. compared with i/^ Vr oscilla- 
tion with maximum mixing. In calculating these 
confidence limits wc have considered correctly the 
non gaussian distribution of the ratio. 

Additional information could be derived from 
the total number of events, at the expense of 
larger theoretical uncertainties. For the best 
value of Ato^ for sterile neutrino oscillation we ex- 
pect a flux reduction of Rfiux = 0.83 for —f Vs 
and Am^ = 0.0025^6^^, to be compared with 
the measured value 0.72. However, due to the 



large theoretical uncertainty, the total number of 
events was not used in the statistical analysis pre- 
sented here. 

It should be noted that this analysis has been 
carried out for the two neutrino mixing case. A 
more complicated oscillation scenario, with 3 or 
more neutrinos fed, or the scenario with large ex- 
tra dimensions [ ^7[ cannot be excluded. 

In conclusion, using the improved statistics af- 
forded by the full MACRO data set, the test of 
the shape of the angular distribution of upgoing 
muons is in good agreement with v^^ — > Vt oscil- 
lation, and maximal mixing. The best is 9.7 
for 9 degrees of freedom. Based on the ratio test, 
the Vs oscillation hypothesis has a 0.033% 

probability of agreeing with the data, and is dis- 
favored at more than 99% C.L.with respect to the 
best fit point of v^^ — > Vt oscillation. 
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